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Abstract 

The information about sizes and nuclear density distributions in unstable (ra- 
dioactive) nuclei is usually extracted from the data on interaction of radioactive 
nuclear beams with a nuclear target. We show that in the case of nucleus-nucleus 
collisions the values of the parameters depend rather strongly on the considered 
theoretical approach and on the assumption about the parametrization of the nu- 
clear density distribution. The obtained values of root -mean-sQuare radii (Rrms) 
for stable nuclei with atomic weights A = 12—40 vary by approximately 0.1 fm 
when calculated in the optical approximation, in the rigid target approximation, 
and using the exact expression of Glauber Theory. We present several examples of 
R-rms radii calculations using these three theoretical approaches and compare these 
results with the data obtained from electron-nucleus scattering. 



PACS numbers: 21.10.Gv, 25.60.Dz, 25.70.Bc 
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1 Introduction 



The root-mean-square nuclear matter radii (Rrms) and the density distributions contain 
an important insight on nuclear potentials and nuclear wave functions. 

In the case of stable nuclei the most precise information comes from the data on 
elastic scattering of fast particles on nuclear targets, ref. [1]. The comparison of the data 
on electron and proton elastic scattering gives separate information about proton and 
neutron distributions in the nucleus, refs. [21 [3]. The Glauber Theory, refs. [HE], [6], is 
typically used to analyze experiments on interactions with nuclei at energies higher than 
several hundred MeV. The obtained values of parameters for nuclear matter density and 
for charge density distributions are summarized in refs. [21[2l[7|. 

Currently, the only way to study density distribution of the unstable (radioactive) 
nuclei is to scatter them off a target. There are a number of technical difficulties in 
the use of hydrogen as target for these experiments, so targets are more often made 
out of different stable nuclei (for example, ^^C, refs. [HI E])- In this case the analytical 
calculation of all Glauber diagrams for nucleus-nucleus interactions becomes impossible 
and some approximative approaches have to be used. The results for nuclear radii and, 
possibly, for nuclear density distributions depend on the considered theoretical approach. 
These problems were already considered for ^He and ^^Li scattering on ^^C target in 
refs. [ini[II]. 

In the present paper we extract the values of nuclear Rrms from the data refs. [SI E] 
on interaction cross sections of stable nuclei with atomic weights A = 12—40 colliding 
with a carbon target. First, we show that the obtained values of Rrms depend on the 
parameterization of the nuclear matter distribution that is used. Then, by assuming 
that the nuclear density in this range of atomic weights can be described by the Woods- 
Saxon expression, we extract the values of Rrms radii from the same experimental data in 
the frame of three different theoretical approaches: the optical approximation, the rigid 
target approximation, and the exact expression of the Glauber Theory. Calculations in 
the Glauber Theory framework are performed using a Monte Carlo simulation technique. 
Finally, we show that the accounting for the non-zero range of NN interaction significantly 
changes the Rrms values. 

Our goal is to demonstrate the numerical level of uncertainties which can appear in 
the analysis of the interaction of unstable nuclei. In general, the optical approximation 
and the rigid target approximation result in smaller values of Rrms when compared with 
the values extracted from the electron scattering data, ref. [7]. The values obtained in 
the framework of the exact Glauber Theory expression are in better agreement with the 
electron scattering data. 
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2 Elastic nucleus-nucleus scattering in the Glauber 
theory 

In the Glauber theory the elastic scattering amplitude of nucleus A on nucleus B with 
momentum transfer q can be written in the frame where -B-nucleus is a fixed target as 
inrefs. [121 [13]: 

FUq)='7^ fd'be^'^'-[l-SABm, (1) 
Zn J 

where k is the incident momentum of one nucleon in A-nucleus in laboratory frame, b is 
the impact parameter, and 



SAB{b) = {A\{B\{l[[l-TMb + u.-s,)]\ \B)\A), (2) 



with 



ih + Ui- s,) = ^Jd'q e-''(^+--^.) . f^^(q) . (3) 

Here Ui and Sj are the transverse coordinates of nucleons, and fNN{(l) is the amplitude 
of elastic nucleon-nucleon scattering, that can be parameterized as 

/^^(9) = ^expf-WV (4) 



47r " V 2 

Here a is the total NN cross section, j3 is the slope parameter of NN elastic scattering. 
We neglect the real part of fNN{(l) since it gives neghgible (~ 1%) contribution to the 
reaction cross section. In the following calculations we use the values a = 43 mb and f3 
= 6 GeV-2. 

Contrary to the case of hadron-nucleus interaction, we can not integrate the Eq. (2) 
anal5d;ically, even with the standard assumption that nuclear densities p(ri,...,r^) in 
both A and B nuclei are the normalized products of one-nucleon densities p(rj): 

A B 

p{ri,...,rA) = Y[pA{ri) , pin, ...,rB) = JlpBiri) , j Srip{ri) = I . (5) 

1=1 i=l •' 

To make the problem manageable, one can retain only part of all contributions in the 
expansion of the product in Eq. (2), that corresponding to the contributions characterized 
by large combinatorial factors. The leading graphs correspond to the so-called optical 
approximation, ref. [2], in which one sums up the contributions with no more than 
one scattering for each nucleon. In other words, only those products of amplitudes 
NN{b + Ui — Sj) in Eq. (2) are taken into accounts which have different indices i, j. This 
approximation corresponds to the summation of diagrams shown in Figs, la, lb, Ic,... 
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In this approximation diagram which describes n-fold interaction has a combinatorial 
factor fyCn ■ bCji- To avoid crowding of hnes in Fig. 1 we have only shown the nucleon 
participants from the nucleus A (upper dots) and B (lower dots) with the links standing 
for interacting amplitudes, and we have not plotted the nucleon-spectators. 



A 











J 



k 



Fig. 1. Diagrams of the interaction of two nuclei, A and B, taken into account 
in different approximations of the multiple scattering theory. 



In the optical approximation averaging and {B\...\B) of the product WieA [1- 

^NN{b + Ui — Sj)] can be substituted with the averaging of the [1 — rjviv(^ + Ui — Sj)]: 



STsib) = 11(^1(^1 [1 - ^NN{b + u, - sj)] \B)\A) . 



Using the standard assumptions of the multiple scattering theory one obtains: 



1 - \T^t{b) 



where 



Topt{b) 



a 



d bid b2TA{bi) ■ Tb(62) • exp 



ex^[-Topt{b)\ , 

(6 + 61-62)' 
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(6) 

(7) 

(8) 



with 

T^(6) = A dzpA {^JWn^) . (9) 

Neglecting the NN interaction range in comparison to the nuclear radii, we have 

Topt{h) = J (fhTAih - hi) ■ TBih) . (10) 

In the diagramatic language Eq. (2) accounts for all possible intermediate states of 
nucleons between the interactions, as it is shown in Fig. 2a, while the optical approxima- 
tion, Eq. (6), would correspond to the interactions with only one pole (nuclear ground 
state) in the both A and B intermediate states, ref. [I3] (see Fig. 2c). 




Fig. 2. Two- fold interaction of two nuclei in the multiple scattering theory (a), in the 
rigid target approximation (b) and in the optical approximation (c). 



Unfortunately, numerical calculations in ref. [15] (see also ref. [16] for the case of 
collisions of very light nuclei) demonstrate that the optical approximation is not accurate 
enough even for the integrated cross sections. The difference with the data amounts for 
~ 10 — 15 % in 0"^^ and it is even greater for differential cross sections, ref. [15]. This 
disagreement can be explained by the fact that series with smaller combinatorial factors 
in Eq. (2) give significant global corrections to the optical approximation results. As a 
matter of fact, the terms of the series are alternating in sign, so, due to the cancelations 
of terms with opposite signs, the final sums of these series can take very different values. 
Thus, some classes of diagrams with non-leading combinatorial factors give significant 
contributions to the final total value. 

The rigid target (or rigid projectile) approximation, refs. [TTl [18], is more explicit 
than the optical approximation. It corresponds to averaging {B\...\B) inside the product 
in Eq. (2): 

SlUb) = {A\ |n - ^NN{b + u,- s,)] \B)\ \A) = [Tr.,Xb)f , (11) 
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where 



Tr.Ab) 



A 



1 



/ 



(fbiTA{bi) exp 



47r/3 



a 



I 



d%2TB{b2) exp 



{b + bi 



2(3 




(12) 



Neglecting the NN interaction range in comparison to the nuclear radii we have: 



This approximation corresponds to the sum of the diagrams in Figs, la, lb, Ic, ... 
and the correction diagrams in Figs. Id, le. If, Ig,... Diagrams in Figs. Id, le. If, Ig,... 
represent the case when each nucleon from the nucleus A can interact several times, 
but all interacting nucleons from B are still different. Each correction diagram which 
describes n-fold interaction has a combinatorial factor smaller than the combinatorial 
factor leading diagrams aCu bCu (Fig- la, lb, Ic, ...). Although, due to the obvious 
asymmetry in contributions of the two nuclei such approach can be theoretically justified 
in the limit A/B -C 1, i.e. for C-Pb, or S-U collisions, this approximation can be used 
sometimes in the case of heavy ion collisions with equal atomic weights. 

Further corrections to the elastic amplitude (some of them are shown in Figs. Ih, li, 
Ij, Ik,...) have been considered in refs. [TU [20l [211 122] • However, the results of such 
corrections are rather complicated for practical use. 

The possibility to obtain the Glauber Theory results without any simplification by 
the direct calculation of Eq. (2) using Monte Carlo simulation was first suggested in 
refs. [231 [21] • This method was used for numerical calculations in refs. [TDl El [25] . 
The simplest algorithm considering the values of coordinates uniformly distributed in 
the interaction region can not be applied here because in most cases several coordinates 
have values corresponding to very small nuclear density. 

The algorithm proposed by Metropolis et al in [26] allows to generate a set of co- 
ordinates which re distributed according to a pre-defined distribution. The Metropolis 
method is as follows: 

1. The initial coordinate Si is randomly generated from the appropriate interval 

2. To obtain next coordinate, a shift Asj is randomly generated and then added to 
the initial coordinate 

3. New coordinate is accepted when the ratio r = p{si + Asj)/p(si) > 1 

4. If the ratio r < 1, new coordinate is accepted only when r > x, where x is a new 
random number x from [0, 1] interval. Otherwise new coordinate is not accepted 

Generated set of coordinates was used to calculate average value of — T^^i^b + 
Ui — Sj)]. Finally, we calculated SAsib) with several different sets of nucleon coordinates. 



Tr.gXb) = jI d%,TAib,-b)-exp [-|Tb(6i) 



(13) 
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3 Monte Carlo simulation of reaction nucleus-nucleus 
cross sections 

The total inelastic (reaction) cross section for the collisions of nuclei A and B, o"^^, is 
equal to the difference of the total interaction cross section cr^^ and the integrated elastic 
scattering cross section, a^^: 

= <B - <B = J d'b[l - \SAB{b)n . (14) 

where a^^B '^ab 

a% = ^ImFUd = 0) = 2/ d%[l - SABib)] , (15) 



al'^ = J d'b[l-SABib)Y , (16) 
Only interaction cross section cr^^ have been measured in refs. [Sj [9]. The difference 

(r) (I) 

between and is that the reaction cross sections include the cross sections of 
all processes except of the elastic scattering AB — AB, whereas the interaction cross 
sections do not include the processes with a target nucleus exitation or with disintegration 
AB AB* {B* 7^ B), hence, < ctab- However, the difference between and 
has been estimated to be less than a few percent for beam energies higher than 
several hundred MeV per nucleon, refs. [271 128] . That allows us to neglect in the present 
paper the difference between cr^^ and cr^^ and compare our calculations of cr^'^g to the 
experimental data from refs. [HI E] on in the same way as it was done in refs. [HI HI] . 

For the numerical calculations it is necessary to use an expression for the nuclear 
matter density distributions in the colliding nuclei. The most detailed information about 
these distributions comes from the data on the differential elastic scattering cross sections 
on nuclear targets. The nucleon density in the light (A< 20) nuclei can be described by 
a harmonic oscillator (HO) density distribution, ref. [H]: 



A/2-2 fr\'^\ ^( 



p^(r) = pi- (l + ^^- (-) ).exp(-- , (17) 



where A is the nucleus size parameter and pi is the normalization constant, while the 
nuclear density distributions in not very light nuclei can be reasonably described by the 
Woods-Saxon expression 

P^^''^ = TT — fr — • 

1 + exp [{r — c)/a) 

Here pi is the normalization constant, c is a parameter measuring the nuclear size, and 
a is related to the diffuseness of the surface, in other words, to the thickness of the 
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nuclear skin. The parameter c shows the value of r at which p(r) decreases by a factor 
2 compared to p(r = 0), p(r = c) = |p(r = 0). The value of a determines the distance r 
(r = 4aln3 ~ 4.4a) at which p(r) decreases from 0.9p(r = 0) to 0.1p(r = 0). 

From the experimental data on cr^^ it is possible to determine only one parameter of 
nuclear matter distribution, let's say Rms: 

Rrms = , = J T^pA{r)d\ . (19) 

It is needed to note that the value of Rrms is rather smaller than the standard nuclear 
radius Ra — 1.2A^^'^ , fm. For example, in the case of uniform nuclear density with radius 

Ra-, Ra = 

The reaction cross section of ^^C-^^C interaction as a function of the Rrms of the ^^C 
nucleus was calculated in different theoretical approximations for various nucleon density 
distributions. The results of the calculations are shown in Figure 3. In these calculations, 
which correspond to energies of 800-1000 MeV per projectile nucleon, the total NN cross 
section value, averaged over pp and pn interactions, a^Jv = 43 mb was used. The nuclear 
matter distribution parameters A, for harmonic oscillator (HO) density in Eq. (17), and 
c, for the Woods-Saxon density in Eq. (18), have been fitted to obtain the required Rrms, 
whereas the parameter a in Eq. (15) has been fixed to the value a = 0.54 fm. 

As the first step we have reproduced the result from ref. [8] for cri2^_i2c the opti- 
cal approximation without NN interaction range, using the harmonic oscillator density. 
Dependence of the reaction cross section, cr|2^_i2(^, on the Rrms is shown by curve 2 of 
Fig. 3. It is in a good agreement with the result from ref. [8], which is shown by marker. 

Curve 1 in Fig. 3 represents the calculations of the dependence of the reaction cross 
section ail^(^_i2(^ on Rrms which was done in the optical approximation with the Woods- 
Saxon density distribution. By comparing curves 1 and 2 of Fig. 3 one can see that 
equal Rrms with different assumptions about the nuclear density distribution result in 
different reaction cross sections. In other words this equivalently means that the same 
experimental reaction (or interaction) cross section with different assumptions about 
nuclear density distribution result in different Rrms, e.g. the assumption of the Woods- 
Saxon density distribution in ^^C nucleus leads to a smaller value of Rrms than the one 
obtained with harmonic oscillator density distribution. 

The corresponding results obtained in the rigid target approximation are shown by 
curves 4 and 3 in Fig. 3. Here again the assumption of Woods-Saxon density distribution 
results in a smaller value of Rrms than the one calculated with the HO density distribution. 
The rigid target approximation contains additional diagrams Figs. Id, le. If, Ig,..., which 
increase the shadow effects. That explains why both curves 3 and 4 lie below curves 1 
and 2. 

The curves 5 and 6 in Fig. 3 show the results of the calculation of all the diagrams of the 
Glauber Theory by using Monte Carlo method and accounting for the finite range of NN 
interaction. Curves 5 and 6 were calculated with the Woods-Saxon density distribution 
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and the HO density respectively. Here, new shadow corrections of the type shown in 
Figs. Ih, li, Ij, Ik,..., appear by comparison to the rigid target approximation. As a 
result, the calculated reaction cross section becomes now smaller at the same value of 

The values of the interaction cross sections 0"{2^_i2^ presented in ref. are 856 ± 9 
mb and 853 ± 6 mb at energies 790 MeV and 950 MeV per nucleon, respectively. The 
older experimental measurement in ref. [22] gives a value (j[i]j_i2Q = 939 ± 49 mb at 
energy 870 MeV per nucleon, i.e. a little larger cross section. The total ^^C-^^C cross 
section was measured to be 1254 ± 54 mb at the same energy , whereas the Glauber 
Theory with Rrms taken from data in ref. |8] predicts a value (j12c_i2c = 1405 mb. 

The Rrms values extracted from the measurement of the interaction cross sections 
of stable projectile nuclei at energies 800—1000 MeV scattered from ^^C target (see 
refs. in [9]) are shown in Table 1. The Rrms values were calculated assuming that the 
nuclear matter density distribution can be described by the Woods-Saxon expression 
Eq. (18). To get the dependence of the reaction cross section a^"^^ on the Rrms, we varied 
parameter c of the density distribution and kept parameter a as a constant at a = 0.54 
fm. The Rrms values were extracted from the agreement of the calculated value a^^^ with 
the experimental values of a^^\ 

In Table 1 one can see that our calculations in the optical approximation with the 
Woods-Saxon density distribution and neglecting the NN interaction range result in the 
shghtly smaller values of Rrms (0.05—0.1 fm) than those obtained in refs. [8], [9]. Rrms 
values are getting even smaller when calculated with a finite range of NN interaction. 
In the case of the Glauber Theory with the Woods-Saxon density distribution shadow 
corrections lead to larger values of Rrms than in the other calculations. 



Nucleus 


Without NN 


range 


With NN range 


Glauber Theory 


HO, optical, [a [9] 


WS, optical 


WS, optical 


WS, rigid target 




2.31 ±0.02 


2.25 ±0.01 


2.09 ±0.01 


2.18 ±0.01 


2.49 ±0.01 




2.47 ±0.03 


2.42 ±0.03 


2.23 ±0.03 


2.35 ±0.04 


2.64 ±0.03 


016 


2.54 ±0.02 


2.48 ±0.02 


2.29 ±0.02 


2.41 ±0.03 


2.69 ±0.02 


pl9 


2.61 ±0.07 


2.55 ±0.08 


2.34 ±0.08 


2.44 ±0.09 


2.75 ±0.07 


Ne^o 


2.87 ±0.03 


2.84 ±0.04 


2.63 ±0.03 


2.75 ±0.04 


2.99 ±0.03 


Na23 


2.83 ±0.03 


2.73 ±0.04 


2.52 ±0.04 


2.62 ±0.04 


2.91 ±0.03 


Mg24 


2.79 ±0.15 


2.65 ±0.23 


2.44 ±0.22 


2.53 ±0.24 


2.85 ±0.20 


C135 


3.045 ±0.037 


2.92 ±0.04 


2.68 ±0.04 


2.76 ±0.04 


3.08 ±0.04 


Ar^o 


3.282 ±0.036 


3.16 ±0.04 


2.90 ±0.03 


2.98 ±0.04 


3.30 ±0.03 



Table 1. The values of Rrms extracted from the measurements refs. [8^ Jj of interaction 
cross section in collisions of projectile nuclear beam with ^^C target at energies of 800- 
1000 MeV per nucleon. 
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4 Comparison of Rrms obtained from nucleus-nucleus 
collisions to electron scattering data 

Let us compare the values for Rrms we have obtained from nucleus-nucleus collisions 
with the published results. It is known from refs. P| [5] that radii of proton and neutron 
distributions in nuclei with Z ~ ^4/2 are practically equal, so we can compare calculated 
radii for nuclear matter with electrical charge radii presented in ref. [7]. 

It is necessary to make a distinction between the distributions of the centres of nu- 
cleons pa(^) and the folded distributions PA(r), where density pAiy^ is convoluted with 
the matter or charge density in the nucleon, Pn{t\- 

PAir) = J pA{r - ri)pM{ri)Sri . (20) 

We have to deal with pa('") and with p^(r) when we calculate Rrms with and without 
accounting for the range of NN interaction, respectively. The r.m.s. radii of pa{t) and 
Pa(?"), Rrms and Rrms, are different and the following relation between them was used in 
ref. 0: 

^Ims = RLs + (0.82 fm)^ . (21) 

In Fig. 4 we compare the result of Rrms calculations in the optical approximation 
(triangles) and in the rigid target approximation (squares) with Rrms values extracted 
from the electron-nucleus scattering experiment, ref. [7]. Calculations of Rrms in both 
approximations were done with zero range of NN interaction (i.e. the folded distribution 
PA^r) was used) and Woods-Saxon density distribution. Obtained results are system- 
atically smaller than the data presented in ref. |7|. This supports our point of view 
that in the optical approximation and in the rigid target approximation the effects of 
nuclear shadowing are too small, and, thus, one obtains agreement to the experimental 
nucleus-nucleus cross section with a smaller value of Rrms • 

The same calculations of Rrms were done in the framework of Glauber Theory. The 
results of these calculations are presented in Fig. 5. In this case it is impossible to provide 
calculations with zero range of NN interaction because the contributions of diagrams with 
loops (see for example. Fig. Ik) present this range in the denominator. 

The Rrms values calculated with distribution pA^r) together with the Rrms values 
extracted from electron-nucleus scattering experiment are presented in the left panel of 
Fig. 5. They are in slightly better agreement to electron-nucleus scattering data than in 
the case shown in Fig. 4. 

However, as we discussed before, the electron scattering data are related to the folded 
distributions pa('")B 

In order to make a more reasonable comparison, we calculate the needed Rrms in 
the case of the Glauber Theory as Rrms + A, where A was calculated as a difference 

^We neglect the difference of electromagnetic and strong interaction nucleon radii. 
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between the Rrms obtained with the distribution of the nucleon centers (see for example 
Eq. 18) and Rrms obtained with the nuclear matter density distribution, Eq. 20. Both 
root-mean-square radii were calculated in the optical approximation. 

The results are shown in the right panel of Fig. 5. The agreement for not very light 
nuclei (A > 20) seems to be reasonable. At A < 20 the Glauber Theory overestimates the 
radii, what probably means that the assumption about Woods-Saxon density distribution 
is not good enough for these nuclei. One can see from Fig. 3 that using of the HO density 
distribution results in smaller r.m.s. values for a given nucleus-nucleus cross section. 

5 Conclusion 

We have shown (see Fig. 3) that the calculations of reaction cross sections for nucleus- 
nucleus collisions with the same Rrms, but different assumptions about nuclear density 
distribution leads to different results. 

These results depend on the used theoretical approximation, the numerical differences 
among the values of reaction cross sections for nucleus- nucleus collisions obtained with 
the same nuclear density distribution in the optical approximation, in the rigid target 
approximation, and in the Glauber Theory, being rather significant, ~ 0.05 — 0.1 fm. 

The nuclear radii obtained from nucleus-nucleus interactions in the Glauber Theory 
are in better agreement with the electron- nucleus data than the ones obtained in the 
optical approximation and in the rigid target approximation. 

We are grateful to G.D. Alkhazov for useful discussions. This work was supported by 
Ministerio Educacion y Ciencia of Spain under project FPA 2005-01963 and by Xunta 
de Gahcia (Spain). It was also supported in part by grant RSGSS-1 124.2003.2 and by 
NSF EPSCoR Program (RSF-023-05). 
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R , fm 

r.m.s.' 



Fig. 3. The cross section of the reaction ^^C— ^^C as a function of the Rrms of the ^^C nu- 
cleus calculated in different approximations with two different nucleon densities: 1) the 
optical approximation without range of NN interaction with Woods-Saxon density distri- 
bution; 2) the optical approximation without range of NN interaction with HO-potential 
density distribution; 3) the rigid target approximation without range of NN interaction 
with Woods-Saxon density distribution; 4) the rigid target approximation without range 
of NN interaction with HO-potential density distribution; 5) Glauber calculation with 
Woods-Saxon density distribution; 6) Glauber calculation with HO-potential density 
distribution. Calculation in ref. [8] is shown as a marker (•) for comparison. 
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Fig. 4. The values of Rrms extracted from electron-nucleus scattering ref. [Tj (filled 
circles), from nucleus-nucleus collisions in the optical approximation (triangles), and in 
the rigid target approximation (squares), both with zero range of NN interaction and 
with Woods-Saxon density distribution. 
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Fig. 5. The values of Rrms extracted from electron- nucleus scattering ref. [7] (filled circles) 
and from nucleus-nucleus collisions in the Glauber Theory (open square) with Woods- 
Saxon density distribution and with finite range of NN interaction (left panel) and with 
zero range of NN interaction (right panel). 



16 



